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Understanding strongly correlated quantum systems is a central problem in many areas of physics. 
f**"» ■ The collective behavior of interacting particles gives rise to diverse fundamental phenomena such as 

confinement in quantum chromodynamics, phase transitions, and electron fractionalization in the 
quantum Hall regime. While such systems typically involve massive particles, optical photons can 
also interact with each other in a nonlinear medium. In practice, however, such interactions are 
■ often very weak. Here we describe a novel technique that allows the creation of a strongly correlated 

' quantum gas of photons using one-dimensional optical systems with tight field confinement and 

, coherent photon trapping techniques. The confinement enables the generation of large, tunable 

optical nonlinearities via the interaction of photons with a nearby cold atomic gas. In its extreme, 
we show that a quantum light field can undergo fermionization in such one-dimensional media, which 
, can be probed via standard photon correlation measurements. 

5_i ' 

qj ' The idea of using nonlinear effects to create optical systems with unusual properties has been pursued for several 
,-£h . decades, and fascinating advances have made in recent years [H, 0, [3, 0, [3] • The focus of such work has mostly been on 
systems containing a small number of photons. At the same time, the effects of strong correlations manifest themselves 
most dramatically in many-body systems, often resulting in new states of matter with properties that are very different 
from those of the underlying particles. One famous example is the Tonks-Girardeau (TG) regime of interacting bosons 
in one dimension, in which strong interactions between particles lead to an effective "fermionization" of bosons [3, 0]- 
This unusual regime has been explored in several condensed matter systems and in recent experiments with ultracold 
atoms in one dimensional traps [3, [3] ■ In this article we investigate the feasibility of creating and detecting a TG 
gas of photons. This system would correspond to nonlinear quantum optics in its extreme, in which individual 
photons behave as impenetrable particles. In this limit, an optical pulse separates into non-overlapping wavepackets 
of individual photons, and a "crystal of photons" can be created. 

Several papers have recently considered the possibility of phase transitions involving photons, focusing on large 
J> ' systems of coupled optical cavities 13, 11 , 13, HI • The present study extends this important work along two principal 
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directions. First, while dramatic progress has been achieved in controllin g in dividual atoms and photons in single 
cavities [3, 13, 0] , the complex architecture of coupled cavities proposed in 13, II , 13, 13] represents a considerable 
experimental challenge. In contrast, the approach described here to realize a TG gas of photons involves currently 
available experimental techniques. Second, we show how signatures of a strongly interacting photon gas can be 
detected using standard quantum optical measurements. This potentially enables many new applications and allows 
' one to study even more exotic phenomena involving interacting particles using readily controlled photons with tunable 
C^*) interactions. 

STRONGLY INTERACTING PHOTONS IN ONE DIMENSION: THE SYSTEM 

Recently, much effort has been directed toward realizing single-mode optical waveguides where the guided photons 
can be tightly confined in the transverse directions to an area A e ff near or below the diffraction limit. Such confinement 
is desirable in part because it allows for a large interaction strength between single photons and nearb y co upled atoms. 
Specific systems that have recently been explored in this context include ta per ed optical fibers [14[, hollow-core 
photonic crystal fibers 13, 13] , and surface plasmons on conducting nanowires |l7t Il8j . 



The propagation and interaction of photons in such a medium can be controlled by interfacing them with atoms 
using quantum optical techniques such as Electromagnetically Induced Transparency (EIT) [f|. In particular, our 
scheme relies on resonantly enhanced optical nonlinearities with low losses [19| using EIT and the creation and 
trapping of stationary pulses of light in the medium using spatially modulated control fields [13] (see Fig. QJi). As 
will be shown, the dynamical evolution of the photonic system is governed by an equation that has the form of the 
Nonlinear Schrodinger Equation (NLSE), where the signs and strengths of the effective mass and interaction can be 
controlled using external fields. Under conditions where the nonlinear interaction is large and effectively repulsive, 
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the TG regime of photons can be achieved. 

Specifically, we consider the propagation of counter-propagating quantum fields, characterized by operators E±, 
inside a single-transverse mode optical waveguide and interacting with cold atoms in the four-level configuration 
shown in Fig. [TJd. The fields E± couple the ground state \a) (in which the system is initialized) to excited state \b) 
with a strength given by g, while metastable state |c) and \b) are coupled by classical, counter-propagating control 
fields fl±(t) (we assume these fields are also guided). Here, the lambda configuration consisting of states \a,b,c) 
comprise the typical EIT setup. In particular, in the case where the quantum and classical fields propagate only in 
one direction, the quantum field can be dynamically and reversibly mapped into a stationary spin-wave excitation by 
turning the control field Q(i) to zero adiabatically. On the other hand, as shown in [2(|, by creating a standing wave 
formation with the control fields (£l+(t) = f2_(i)), the quantum fields can be effectively stopped while simultaneously 
maintaining a non-zero photonic component of the excitation (see Fig. QJl) . Intuitively, the standing wave pattern 
forms a set of reflection gratings that trap the photonic excitations through multiple scattering. The nonlinear response 
is introduced via an additional state \d) that is coupled to |c) by the quantum fields (for simplicity we assume that 
this coupling strength is also given by g). When these fields are far off resonance, this coupling results in an ac Stark 
shift of state |c) whose magnitude is proportional to the intensity of the quantum field 19]. This gives rise to an 
intensity-dependent refractive index, or nonlinear susceptibility, that in turn influences the evolution of the quantum 
fields. While similar schemes for realizing nonlinear optics in atomic media have been previously explored [211 ] . their 
implementation in waveguides with tight confinement is unique because they can constitute a true one-dimensional 
system, and because the strong transverse localization enables large nonlinear interactions. 



THE LIEB-LINIGER MODEL WITH STATIONARY PULSES OF LIGHT 



Considerable theoretical literature exists discussing the connection between Maxwell's equations for a nonlinear 



medium in one dimension and the Lieb-Liniger Model (LLM) describing interacting massive particles 2 2l. 1 2 31] . However, 
most of this work has focused on the regime of attractive interactions. Interest in this regime stems from the 
observation that interaction with basic two-level atomic systems leads to photon bunching, since it is easier for many 
photons to pass through an atom without being absorbed [24| . The main feature of the attractive regime is the 
formation of solitons [22j. This regime of the LLM has few quantum mechanical features and can be described 
by classical equations of motion. Indeed, typical optical solitons contain on the order of a million photons, and 
therefore a quasi-classical description is very appropriate [251 ]. On the other hand, the repulsive case is intrinsically 
"quantum mechanical" , as one needs to keep track of correlations at the level of individual particles. Solutions based 
on perturbation theory or classical equations of motion break down and one needs to use non-perturbative approaches 



such as the Bethe ansatz solution [261 1271. 1281. 1291] or the Luttinger liquid formalism [3fJ. Thus, formation of a TG gas 
of photons is fundamentally a collective many-body effect. 

We now derive an evolution equation for the fields of the system illustrated in FigureQ] Following the methods of [20L 
we define dark-state polariton operators 'f ± describing the coupled photonic and spin- wave excitations, which 
in the slow-light limit are given approximately by — g^/2 / nn z E±/Vt± 1 where n z is the density of atoms coupled 
to the waveguide (the density is assumed to be uniform). We specialize to the case when f2±(t) = fl(t). We further 
assume that the quantum and control fields vary slowly in time, such that the fast-varying atomic operators can be 
adiabatically eliminated, while the remaining slowly-varying operators are solved in the adiabatic limit. Subsequently 
inserting these solutions into the Maxwell-Bloch equations describing evolution of the quantum fields yields an effective 
NLSE for the photon gas (see Methods), 

i8t9(z,t) = -^B%*(z,t) + 2gtf(z,t)* 2 (z,t), (1) 

where 

+ n*—g±, 2 S= r -lp. (2) 

2 4A w 9 A p 

Here Tid = 4:irg 2 /v is the spontaneous emission rate of a single atom into the waveguide modes, where v is the 
velocity of these modes at the atomic resonance frequency in an empty waveguide, while v g ~vQ 2 / (irg 2 n z ) is the 
group velocity of untrapped pulses under EIT conditions. Ao and A p are the one-photon detunings of the fields E± 
from the transitions |o)-|6) and |c)-|d), respectively (see Fig. []})), and <7~l/\A4cff is the single-photon, single-dipole 
interaction matrix element. In principle the full dynamics of the field will also include losses and higher-order terms, 
and the conditions under which such terms can be neglected are described in Methods. 
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Eq. (p} determines the evolution of a quantum field $f(z, t) as derived from the Hamiltonian of the Lieb-Liniger 
model, 

H = h [ dz\-^— d z ^(z)d z t>(z) + g^(z)^(zW(zW(z)}. (3) 
J 2m cff 

The first term on the right describes the kinetic energy in one dimension of bosons with mass m e s, while the second 
term describes a contact interaction potential. The quantum field *B(z,t) satisfies the usual equal-time bosonic 
commutation relations, [^(z, t), &<(z', t)] = S(z — z'). The behavior and properties of this system can be effectively 
characterized by a single dimensionless parameter 

-■ = m ° ff g = n * , (4) 

n ph 8A A p n ph ' 

which physically corresponds to the ratio of the interaction and kinetic energies. Here n p h is the density of photons at 
the center of the pulse. When 7 < 0, the interaction with the atoms induces an effective attraction between photons. 
As discussed above, this is responsible, e.g., for the formation of bound states (solitons) when a large number of 
photons are present. On the other hand, for 7 > the regime of effective repulsion between photons is realized. The 
special limit 7 — * 00 is called the Tonks-Girardeau (TG) regime 

0,0- From E q- ©> 

one finds that the regime of 

repulsion can be achieved if exactly one of the detuning parameters Ao or A p is negative. It is also clear that only 
the attractive regime can be reached using a two-level system, since Ao = A p in this situation. Finally, we note 
that jcxTfj-) oc A~g , which underscores the importance of tight mode confinement in reaching the strongly interacting 
regime. 

One useful feature of this realization of the NLSE is that the parameters m e e and g can be dynamically tuned by 
varying different parameters of the system. Specifically, both are functions of detuning, and thus 7 can be altered 
by changing either the control field frequencies or by externally manipulating the energies of the atomic levels. As 
described below, this tunability facilitates the creation of novel photonic states. 



PREPARATION AND DETECTION OF STRONGLY CORRELATED PHOTON GAS 



The process by which novel photonic states can be prepared and detected consists of three basic steps - loading 
of the pulse, controlled evolution under the NLSE, and readout of the final photonic state - each of which we now 
describe. During the initial loading process, a resonant optical pulse, given by, e.g., a coherent state, is incident from 
one direction. It is injected into the waveguide at the same time that the co-propagating control field (say f2 + (£)) is 
turned on. During the loading procedure the counter-propagating control field f2_(f) is off. This describes the usual 
situation in EIT, where the input field becomes spatially compressed upon entering the medium and propagates with a 
variable group velocity v g ~vQ.+ (t) / (g 2 n z ). Once the pulse completely enters the medium, f2+(i) is adiabatically turned 
to zero, reversibly converting the photonic excitation into a pure spin- wave excitation in the atomic medium [3lj ]. 
Under certain conditions this input pulse can be stored with minimal distortion relative to the initial field [31j . such 
that all relevant properties (e.g., field correlations) remain constant during the loading. Following the initial storage, 
both control fields are then adiabatically switched on, with Q± — O(t), during which the pulse becomes trapped and 
evolves under Eq. |T]). The parameters m e a and g can be changed in time during the evolution to reach the final 
state of interest. Once this state has been achieved, the pulse is released by turning one of the control fields (say 
f2_) off, thereby allowing the pulse to propagate undistorted until it exits the waveguide During this readout 
of the pulse, any spatial correlations that formed while evolving under the NLSE are directly mapped into temporal 
correlations (at a common point in space) in the outgoing field, which can be measured using standard quantum 
optical techniques. 

A characteristic signature of a strongly interacting gas is the appearance of Friedel oscillations [32] in the normalized 
second-order correlation function, g^ 2 \z,z') = (I(z)I(z'))/((I(z))(I(z'))), where I(z) — ^ (z)^ (z) is the stationary 
pulse intensity prior to release. In particular, if the photonic state is close to the ground state of the LLM, g^ 2 \z, z') 
contains an oscillating part that behaves as ~ cos(2kp(z — z')) (see Methods), where kp = irn p h is the Fermi 
momentum in the TG limit. In the TG limit, the ground-state correlation function takes on a simple form given by 
Lenard's formula [33} , 

^(-')--(l^f) 2 - w 
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The 2feF-oscillations are a direct manifestation of "fermionization" of bosons. We note that Friedel oscillations are 
more than simple anti-bunching in that they indicate real crystal correlations. In particular, one cannot predict 
the position of an individual photon, but knowing the position of one photon, other photons are likely to follow at 
well-defined distances determined by the average photon density. These correlations are predicted to decay relatively 
slowly in space. 

The usual definition of the TG regime of the LL model is given with respect to the equilibrium state. However, in 
our system the evolution of the initial coherent state under the NLSE inherently involves non-equilibrium quantum 
dynamics, and hence one must specify the conditions under which one can obtain a strongly interacting state of 
photons that is close to the ground state of the LLM and exhibits its characteristic features, such as Friedel oscillations. 
Specifically, we consider the evolution of a pulse under Eq. ([I]), where the atomic parameters are varied such that the 
effective mass m e g is constant in time, while the interaction strength g(t) increases exponentially. From Eq. Q, the 
resulting time-dependent interaction parameter can be written in the form 7 = 7oe /3 " F *, where (3 is a dimensionless 
parameter characterizing the rate of increase, and LL>F^n 2 h /m c g corresponds to the "Fermi energy" at the center of 
the pulse before expansion. We assume that 70^1, so that the system initially consists of non-interacting photons, 
and interactions are gradually switched on to reach the regime with 7 ^> 1. Here we will not discuss the effect of a 
time-dependent mass, although it can be included using a similar analysis. We also note that the related problem 



involving the sudden switching of the interaction strength from zero to a large value has been studied in 34 1. 

The non-equilibrium dynamics of the LLM model has previously been studied in the context of ultracold atoms 
and has focused either on changing the interaction strength in a system with uniform, constant density HH, 36| , or on 



the expansion of particles in a system with constant interaction [371 . 1 3 81 ] . In our system both processes take place, but 
in the experimentally relevant regime of a large photon number there is a separation of time scales that simplifies the 
analysis. We consider an initial pulse containing N p h~n p hZo photons with spatial extent zq at t = 0. Our discussion 
applies for a general pulse shape and thus we need not specify it. 

There are three distinct regimes in the time evolution of the pulse: i) Interactions are weak and the photons expand 
freely due to dispersion; ii) Interactions begin to dominate over the kinetic energy. However, the system is still in the 
weakly interacting regime with 7 < 1, in which case the expansion is hydrodynamic M; iii) The system reaches 



the strongly interacting regime with 7 > 1, and the expansion resembles that of fermionized bosons [37j, l38(. A simple 
analysis (sec Methods) shows that by the time the system reaches regime iii), the relative change in the photon density 
at the center of the pulse is only of the order of „ 2 ^ . Hence for (3N p h > > 1 one can assume that the turning on of 

interactions takes place at a constant density and any subsequent expansion takes place in the TG regime. 

We now consider the effect of non-adiabaticity on correlation functions such as g^ 2 \z, z'). By analogy with ultracold 
atoms, we introduce a chemical potential, which in regimes i) and ii) is given by /j,(t) ss g(t)n p h- We can approximately 
separate the turning on of interactions into two stages, which both take place during the regimes i) or ii). In the first 
stage, fi > fj, 2 and the evolution is diabatic. In the second stage, fi < n 2 and the evolution is essentially adiabatic. 
At the time t ac iiab separating the two regimes, the interaction parameter is given by 7(i £2ia&)~/3. The first stage can 
be thought of as an instantaneous projection of the wavefunction, which gives rise to a finite density of excitations 
characterized by an effective healing length £neij~(MW m eff) _1 ^ 2 at time t a di a b- During the second stage, the number 
of excitations does not change. Hence, we find that the finite rate of change in the interaction strength leads to a finite 
correlation length in our system, £ neq ~ P~ 1 ^ 2 ' n p h - F° r length scales shorter than £n eg , ah correlation functions are 
essentially the same as in the ground state, whereas for length scales longer than £n ei j, correlation functions rapidly 
decay. To observe Friedel oscillations over length scales on the order of the inter-photon distance, for example, requires 
that [3 < 1 . The argument presented above can be turned into a quantitative calculation for the correlation functions 
following the time dependent change in the interaction strength [36j |, This analysis uses bosonization to treat the 
LLM and the conclusions agree with the qualitative picture presented here. 

At the end of regime ii) we have a system of "fermionized" hard core photons that should exhibit Friedel oscillations. 
In regime iii), the pulse of "hard-core photons" expands, but such spreading does not lead to the suppression of the 
Friedel oscillations [37], [H| . The problem of expansion of hard-core photons starting from a general pulse shape has 
to be analyzed numerically However, an explicit analytic solution is available for the case of a parabolic pulse shape. 
Under realistic conditions this solution, which is discussed in detail in the Methods section, yields the correlation 
function shown in Fig. [21 

In principle, during the evolution one must also consider the effects of photon losses, which set a maximum evolution 
time t m ax and interaction parameter ^max that can be achieved before a substantial fraction of the initial pulse is 
dissipated. As derived in the Supplementary Information (SI), 

. / //3|A |\ T OD\ 

W-nnn^a^— j 'tfj^Jfr) , (6) 
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where T is the total spontaneous emission rate of states \b) and \d) (assumed to be equal for simplicity), which includes 
the emission rate into the waveguide modes (Tid) as well as emission into non-guided modes (e.g., into free space). 
We have also defined a single-atom cooperativity 77 = Tm/T (?7<1) which describes the rate of emission into the 
waveguide compared to the total emission, and defined the optical depth of the medium, OD = rjZQn z . We note that 
Imax can be improved by increasing the optical depth or cooperativity of the system, and optimized by adjusting the 
detuning |Ao|. The optimal values of j ma x, as functions of optical depth and cooperativity, are plotted in Fig. O for 
parameters (3 = 1, 70 = 0.1, and N p h = 10. One sees that with realistic values of OL>~2000 and ?7~0.2, for example, 
an interaction parameter of 7maa;~10 is possible. While photon losses limit the maximum evolution time, somewhat 
surprisingly, the nonlinear losses may also help in bringing the system closer to the ground state of the LLM model. 
Specifically, these losses predominantly remove states in which two photons are close to each other, which correspond 
to high energy states of the LLM. Deep in the TG regime we can estimate the time it takes for the system to lose the 
high energy states, « gn p h- On the other hand, the photon loss rate is ~ 3- wjr. Thus t en /ti os w 7 _1 and 
for large 7 there is a sufficient time window for the high-energy states to decay before too many photons are lost. 



OUTLOOK 



The above analysis indicates that strongly correlated states of photons can be controllably prepared and observed 
in one-dimensional waveguides. These techniques are made possible through strong coupling between the photons 
and nearby atoms and the use of quantum optical techniques such as EIT, which allow the system to be widely 
tunable. Such photonic states should find numerous applications in various areas of physics. Crystal correlations that 
arise in the fermionized state make it a promising candidate for applications in metrology and quantum information. 
In particular, TG states feature strongly suppressed photon number fluctuations within a given detection interval. 



Such states therefore could be used as an input for sub-shot noise interferometers [41|, |42| , or in extension to schemes 
for quantum computing [43] or quantum cryptography [lij that rely on single photons. Another exciting direction is 
quantum simulation of matter Hamiltonians using optical systems. In the discussion so far, we have considered photons 
with only one polarization. Including photons of different polarizations should be equivalent to adding a spin degree of 
freedom to effective matter Hamiltonians. This opens up exciting prospects for exploring spin charge separation (see, 
e.g., [HI) and modelling exotic spin systems. It is also interesting to note that the level structure of the atoms 
comprising the medium can vary considerably. This fact can be used for the study of strongly-correlated systems with 
non-Abelian symmetry, similar to the ones realized in multi-channel Kondo models and quantum chromodynamics. 
Possible phases and phase transitions in these models are difficult, if possible at all, to realize in matter systems. 
Furthermore, it would be interesting to consider a situation where interactions between the cold atoms create non- 
trivial correlations [46| and their effect on the resulting photonic states. Using light to simulate matter Hamiltonians 
will give a new meaning to the old idea of the particle-wave duality. 

METHODS 
Derivation of NLSE for photons 

The Hamiltonian corresponding to the system described in Fig. [T] is given in a rotating frame by 

H = -hn z J dz[A a bb + A p a dd + gV2^({a ba + a dc )(E+e lkaZ + E_e~ ikoZ ) + h.c.) 

+ {{n + {t)e lk ^ + n-(t)e- lk < z )<j bc + h.c.)}, (7) 

where E± = E±(z,t) are slowly- varying operators describing the quantum fields, and <7ij = crij{z, t) are collective, 
continuous operators describing the average of over atoms in a small but macroscopic region around z. For 

simplicity, we have assumed equal transition matrix elements g ba = g dc = g between the quantum fields E± and the 
transitions |a)-|6) and |c)-|d), where gij~(i\r\j)-\/u)ij/TieoA c ff. Note that we have also included a one-photon detuning 
Aq between the quantum fields and transition |a)-|6) (while maintaining two- photon resonance), whose purpose is to 
provide a group velocity dispersion or "effective mass" term in the field evolution equations. Defining slowly varying 
operators a ab = a ab ,+e lk ° z + a ab ^e^ lk ° z and a cd = a cd ,+e lk ° z + a cd _-e~ lk ° z , the Maxwcll-Bloch equations describing 
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evolution of the fields under H are given by 



! 9 ± j-^ E±(z,t) = ^ l9Uz (a abi ±(z,t) + a cd ,±(z,t)), (8) 



v dt dz 

while the usual Langevin-Bloch equations describing evolution of can be derived following the methods of 2(| 21 



31]. Following these references, we define polariton operators ^>± to describe the collective excitations of field and spin- 



wave coherence a ac that result from coupling with the control fields, which in the relevant limit that the excitations 
are mostly in spin-wave form are given by \&± = g\/2irn z E±/fl±. To proceed further, we adiabatically eliminate the 
Langevin-Bloch equations for the fast-decaying atomic operators (e.g., a a b and <7 ct j), and slowly-varying operators are 
solved in terms of \&±, discarding higher time derivatives in the slowly- varying limit. Plugging these results back into 
Eq. ([8]), and specializing to the case where fl±(t) = il, we obtain evolution equations for the polaritons alone, 

-d t V + d z A = -—d t ^> - ^ (2^ f ^ + A^A) vE + noise, (9) 

v v g v(2A p + il ) v ' 

-d t A + 8 Z V = - y 2 A- y &*A + noise. 10 

v v(i — 2iA ) v(2A p + il ) 

Here we have defined the symmetric and anti-symmetric combinations "J = (W + + , 3 / -)/2 and A = (W + — 4'_)/2, 
and a group velocity v g ~vfl 2 / (irg 2 n z ) under which the pulses would propagate were they not trapped. The total 
spontaneous emission rates T (which include decay into channels other than the guided fiber modes) from states 
\b), \d) are assumed to be identical for simplicity. We note that there are also noise operators, which are associated 
with the dissipative terms in the equations above. Because we are primarily interested in the regime where losses 
are not significant, the specific form of these operators is not important here. With sufficient optical depth, A can 
be adiabatically eliminated, Am(2iAo — T)v(d z ^)/ (4:Trg 2 n z ), where we have assumed that the nonlinear contribution 
~4 r ^A is small. Physically, this result corresponds to a pulse matching phenomenon (47j between the quantum and 
control fields, whereby any imbalance between ^ + and rapidly goes to zero. Substituting the expression for A 
back into Eq. ©, and considering the relevant case where u 9 <t) yields the NLSE with complex effective mass and 
two-body interaction strength, 



(2A + iT)v g v 2 ^g 2 v g f 2 

Wt Ang 2 n z 2 + v(A p + iT/2) 



(11) 



Identifying T\d = 4:irg 2 /v and ignoring the loss terms reproduces the ideal NLSE given in Eq. (fT]). 

We now consider carefully the limits under which Eq. (fT]) well approximates the complete dynamics of the field. 
First, one requires that the ac Stark shift of |c) due to the nonlinear interaction fits within the frequency range where 
EIT is efficient (i.e., within the transparency window), which is conservatively satisfied when n p h/n z <ti\A p \/\T — 2iAo\. 
In addition, requiring that higher-order derivatives of the field be negligible compared to those appearing in Eq. (TTJ 
places a restriction on the maximum wavevector k max of the spin-wave excitation. In the TG regime, for example, 
kmax^nph and one consequently finds that n p h/n z -^.Tiu/\2Ao + iT\. Finally, as discussed further in SI, one must 
also ensure that the loss terms in Eq. (jlip do not cause the dissipation of too many photons, which sets a maximum 
allowed evolution time t max for the system. 



Density-density correlations in one-dimensional system of bosons 



Evaluation of the density-density correlation function in the ground state is challenging even though the Hamiltonian 
in Eq. ([3]) is exactly solvable. However, in the regime of interest where 7 is large, one can obtain an analytic 
expression (2?| (we assume that translational invariance is present and therefore g^> (z\, t\\ z 2 , t 2 ) — g^ 2 \z\ ~ z 2 — 
Az, ti - 1 2 = At)), 



g<~ 2 \Az,At) = l 



K 



4tt 2 t 



ph J \qi\>Trn p h J \qi\<Ttn v h 



Mii-92) cos[{ qi - q 2 )zVK][l + 



ng 



irn ph 



(- 



<?3 - qi 93 - 92 



where K = 1 + 4/7. Note that g( 2 ' decays as t~ x for a range of t > 0. 

For arbitrary interaction strength one can use the exact solution to numerically evaluate g^> 27, 28, 29|. On the 
other hand this numerical solution [28|, [29[ as well as a number of other arguments suggests that an effective 
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description of a one-dimensional gas of bosons, the Luttinger liquid theory, provides very good long-distance, long-time 
behavior. Using this theory, one can demonstrate that the density-density correlation function in the ground state 
exhibits a power-law decay of correlations and 2kp Friedel oscillations, 

o^lAz At) = 1 + K (AZ)2 = {v ° At)2 + Bco < 2k ^ z ) (13) 



where B is a non-universal constant. The interaction parameter K, or Luttinger parameter can be numerically 

extracted from the exact solution. In the limit of strong interactions, K = 1 + 4/7 and tends to 1 in the TG regime. 
Here kp — Trn p h is an emergent Fermi momentum. 

In the case of sudden switch on of the interaction, one can use the completeness of the Bethe-Ansatz wave functions 
basis and expand the initial state over this basis. Since the matrix elements of the density operator in the Bethe 
states are known (49| , computation of the density-density correlation function starting from an arbitrary initial state 



can in principle be performed 34j |. 



Expansion of the pulse 

We consider the expansion of the optical pulse during stages i) and ii) introduced in the Section "Preparation 
and detection of strongly correlated photon gas". In regime i), the interaction energy of particles can be neglected in 
comparison with the kinetic energy as long as g(t)n p h is smaller than the kinetic energy l/m e gZQ. The time t\ at which 
these energies become comparable satisfies toe — log (./V 2 ^ 70), and thus if one chooses parameters such that N^yo > 1 
free expansion can be ignored (for large photon number one can simultaneously satisfy 70^1 such that the system is 
initially non- interacting) . Reg ime ii), which is valid until ^(t = t c ) = 1, can be analyzed using the usual hydrodynamic 
equations (see, e.g., Refs. |39l lioj]). From the equation of motion for the flow velocity, we find v(t) ~ ^^ gft), and 
from the conservation of particle number, n(t) ~ -j^-vif). Assuming that the change in density at the center of the 
pulse is small, these equations can be integrated to yield a relative change of density An vh / 'n v h^{ 1 2 N^ h . 



Next we consider the expansion in the strongly interacting regime. Recent work [37.I38H demonstrated that spreading 
of a pulse of hard core bosons enhances the "fermionic" character of the wavefunction [37| . While generally this problem 
requires numerical simulations, the solution turns out to be extremely simple for an initial pulse density of the form 
n p h{z) = no(l — z 2 /zq) 1 / 2 , which corresponds to hard core bosons released from a parabolic potential [38j |. In this 
case the structure of the density correlations is preserved and there is only a rescaling of length scales [IH • An explicit 
calculation of </ 2 ' for this system is shown in Fig. [2] Another signature of fermionization of an expanding pulse of hard 
core bosons can be observed in the momentum distribution function n(k) (density at wavevector k), corresponding to 
the Fourier transform of the first order correlation function g^(z, z') — (^(z)^ (z 1 )) . After sufficiently long expansion 



it approaches a Fermi distribution [37J, [38| . We emphasize that the qualitative features of our results remain valid for 
various shapes of the initial pulse, including Gaussian and parabolic pulses, although numerical analysis is needed to 
obtain precise answers. 



* Electronic address: demler@physics.harvard.edu 
[1] Boyd, R.W. Nonlinear Optics (Academic, New York, 1992). 

[2] Haroche, S. & Raimond, J.M. Exploring the quantum: atoms, cavities and photons (Oxford University Press, New York, 
2006). 

[3] Miller, R. et al. Trapped atoms in cavity QED: coupling quantized light and matter. J. Phys. B: At. Mol. Opt. Phys. 38, 
S551-S565 (2005). 

[4] Walther, H., Varcoe, B.T.H., Englert, B.-G. & Becker, T. Cavity quantum electrodynamics. Rep. Prog. Phys. 69, 1325-1382 
(2006). 

[5] Fleischhauer, M., Imamoglu, A. & Marangos, J. P. Electromagnetically induced transparency: Optics in coherent media. 

Rev. Mod. Phys. 77, 633-675 (2005). 
[6] Tonks, L. The Complete Equation of State of One, Two and Three-Dimensional Gases of Hard Elastic Spheres. Phys. Rev. 

50, 955-963 (1936). 

[7] Girardeau, M. Relationship between Systems of Impenetrable Bosons and Fermions in One Dimension. J. Math. Phys. 1, 
516-523 (1960). 

[8] Kinoshita, T., Wenger, T. & Weiss, D.S. Observation of a one-dimensional Tonks- Girardeau gas. Science 305, 1125-1128 
(2004). 



8 



[9] Paredes, B. et al. Tonks-Girardeau gas of ultracold atoms in an optical lattice. Nature 429, 277-281 (2004). 
[10] Hartmann, M.J., Brandao, F. G. S. L. & Plenio, M.B. Strongly interacting polaritons in coupled arrays of cavities. Nature 
Physics 2, 849-855 (2006). 

[11] Greentree, A.D., Tahan, C, Cole, J.H. & Hollenberg, L.C.L. Quantum phase transitions of light. Nature Physics 2, 856-861 
(2006). 

[12] Angelakis, D.G., Santos, M.F. & Bose, S. Photon blockade induced Mott transitions and XY spin models in coupled cavity 

arrays. arXiv:quant-ph/0606159 
[13] Neil Na, Y.C., Utsunomiya, S., Tian, L. & Yamamoto, Y. Strongly Correlated Photons in a Two-Dimensional Array of 

Photonic Crystal Microcavities. [arXiv:quan t-ph/0703219 
[14] Nayak, K.P. et al. Optical nanofiber as an efficient tool for manipulating and probing atomic fluorescence. Opt. Express 

15, 5431-5438 (2007). 
[15] Knight, J.C. Photonic crystal fibers. Nature 424, 847-851 (2003). 

[16] Ghosh, S., Sharping, J.E., Ouzounov, D.G. & Gaeta, A.L. Resonant Optical Interactions with Molecules Confined in 

Photonic Band-Gap Fibers. Phys. Rev. Lett. 94, 093902 (2005). 
[17] Chang, D.E., S0rensen, A.S., Hemmer, P.R. & Lukin, M.D. Quantum optics with surface plasmons. Phys. Rev. Lett. 97, 

053002 (2006). 

[18] Akimov, A.V. et al. Generation of single optical plasmons in metallic nanowires coupled to quantum dots. Accepted by 
Nature (2007). 

[19] Schmidt, H. & Imamoglu, A. Giant Kerr nonlinearities obtained by electromagnetically induced transparency. Opt. Lett. 
21, 1936-1938 (1996). 

[20] Bajcsy, M., Zibrov, A.S. & Lukin, M.D. Stationary pulses of light in an atomic medium. Nature 426, 638-641 (2003). 
[21] Bajcsy, M., Andre, A., Zibrov, A.S. & Lukin, M.D. Nonlinear Optics with Stationary Pulses of Light. Phys. Rev. Lett. 94, 
063902 (2005). 

[22] Lai, Y. & Haus, HA. Quantum theory of solitons in optical fibers. II. Exact Solution, Phys. Rev. A 40, 854-866 (1989). 
[23] P. Drummond, Quantum Theory of Nonlinear Optics, [http: / / ww w.physics.uq.edu.au/people / drummond/NOTES /lecture. pdf 
(2001). 

[24] Chang, D.E., S0rensen, A.S., Demler, E.A., & Lukin, M.D. A single-photon transistor using nanoscale surface plasmons. 

Nature Phys. 3, 807-812 (2007). 
[25] Mazets, I.E. & Kurizki, G. How different are multiatom quantum solitons from mean-field solitons? Europhys. Lett. 76, 

196-202 (2006). 

[26] Lieb, E.H. & Liniger, W. Exact Analysis of an Interacting Bose Gas. I. The General Solution and the Ground State. Phys. 
Rev. 130, 1605-1616 (1963). 

[27] Korepin, V.E., Bogoliubov, N.M. & Izergin, A.G. Quantum Lnverse Scattering Method and Correlation Functions (Cam- 
bridge U.P., Cambridge, UK, 1993). 

[28] Caux, J.-S., & Calabrese, P. Dynamical density-density correlations in the one-dimensional Bose gas. Phys. Rev. A 74, 
031605 (2006). 

[29] Caux, J.-S., Calabrese, P. & Slavnov, N.A. One-particle dynamical correlations in the one-dimensional Bose gas. J. Stat. 
Mech., P01008 (2007); 

[30] Haldane, F.D.M. Effective Harmonic-Fluid Approach to Low-Energy Properties of One-Dimensional Quantum Fluids. 

Phys. Rev. Lett. 47, 1840-1843 (1981). Erratum: ibid 48, 569 (1982).' 
[31] Fleischhauer, M. & Lukin, M.D. Dark-State Polaritons in Electromagnetically Induced Transparency. Phys. Rev. Lett. 84, 

5094-5097 (2000). 
[32] Friedel, J. Metallic alloys. Nuovo Cimento Suppl. 7, 287 (1958). 

[33] Lenard, A. One-Dimensional Impenetrable Bosons in Thermal Equilibrium. J. Math. Phys. 7, 1268-1272 (1966). 
[34] Rostunov, T., Gritsev, V. & Demler, E. In preparation. 

[35] Cazalilla, M.A. Effect of Suddenly Turning on Interactions in the Luttinger Model. Phys. Rev. Lett. 97, 156403 (2006). 
[36] Polkovnikov, A. & Gritsev, V. Breakdown of the adiabatic limit in low dimensional gapless systems. arxiv:cond- 
mat/0706.0212. 

[37] Rigol, M. & Muramatsu, A. Free expansion of impenetrable bosons on one-dimensional optical lattices. Mod. Phys. Lett. 
B 19, 861-881 (2005). 

[38] Minguzzi, A. & Gangardt, D.M. Exact Coherent States of a Harmonically Confined Tonks-Girardeau Gas. Phys. Rev. Lett. 
94, 240404 (2005). 

[39] Castin, Y. Dum, R. Bose-Einstein Condensates in Time Dependent Traps. Phys. Rev. Lett. 77, 5315-5319 (1996). 

[40] Kagan, Yu., Surkov, E.L. & Shlyapnikov, G.V. Evolution of a Bose-condensed gas under variations of the confining 

potential. Phys. Rev. A 54, R1753-R1756 (1996). 
[41] Dowling, J. P. Correlated input-port, matter-wave interferometer: Quantum-noise limits to the atom-laser gyroscope. Phys. 

Rev. A 57, 4736-4746 (1998). 

[42] Bouyer, P. & Kasevich, M.A. Heisenberg-limited spectroscopy with degenerate Bose-Einstein gases. Phys. Rev. A 56 , 
R1083-R1086 (2002). 

[43] Kok, P. et al. Linear optical quantum computing with photonic qubits. Rev. Mod. Phys. 79, 135-174 (2007). 
[44] Gisin, N., Ribordy, G., Tittel, W. & Zbinden, H. Quantum cryptography. Rev. Mod. Phys. 74, 145-195 (2002). 
[45] Schulz, H.J. Fermi liquids and non-Fermi liquids. In Mesoscopic Quantum Physics, Proceedings of Les Houches Summer 
School LXI (Elsevier, Amsterdam, 1995), p. 533-603. Eds: Akkermans, E., Montambaux, G., Pichard, J. & Zinn- Justin, J. 
[46] Masalas, M. & Fleischhauer, M. Scattering of dark-state polaritons in optical lattices and quantum phase gates for photons. 



9 



Phys. Rev. A 69, 061801(R) (2004). 
[47] Harris, S.E. Electromagnetically induced transparency with matched pulses. Phys. Rev. Lett. 70, 552-555 (1993). 
[48] Giamarchi, T. Quantum Physics in One Dimension (Oxford UP, New York, 2004). 

[49] Slavnov, N.A. Nonequal-time current correlation function in a one-dimensional Bose gas. Theor. Math. Phys. 82, 273-282 
(1990). 

We gratefully acknowledge support from the NSF, Harvard-MIT CUA, DARPA, Air Force, and Packard Foundation. 
Correspondence and requests for materials should be addressed to E.A. Dcmler. 



10 




Q + (t) Q_(t) 



FIG. 1: Illustration of fields and atoms comprising the system a) Schematic of fields inside the waveguide, whose axis 
of propagation corresponds to the horizontal axis. The control beams fi±(t) (shown in blue) create a standing wave inside the 
waveguide, which forms a Bragg grating that traps a quantum optical field inside the medium (intensity (E'E) shown in red). 
The optical field couples to spin-wave excitations in the medium, resulting in collective polariton excitations whose density 
is also plotted, b) Schematic of the four-level atomic configuration and coupling between levels and fields used in our 
system. 




FIG. 2: Density-density correlation function g^ 2 \z,z' = 0) for an expanding Tonks-Girardeau gas of photons with initial 
density profile n p h(z) = no(l — z 2 /zq) 1 ^ 2 . This expansion is equivalent to the problem of a TG gas released from an initial 
parabolic confining potential, z' = denotes the center of the pulse, and distances are indicated in units of kp 1 . The density- 
density correlation function shown here is for a system of N p n = 10 photons, zo^5kp , and at a time t — Wujp following the 
initial release. 
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FIG. 3: Maximum interaction parameter "/ ma x as functions of optical depth and single-atom cooperativity, optimized over 
the detuning Aq. 7 mra is plotted for fixed values of /3 = 1, 70 = 0.1, and N p h = 10. It is evident that the TG regime for 
photons can be approached by increasing either the cooperativity or optical depth. 



